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1 Introduction 



The computation of correlation functions in field theoretical models is a difficult problem. 
In this article we present a novel approach which applies to models where a two point 
function can be written as 

S(x,y) = j \p + Q]-J dpt(Q) . (1.1) 

Here P is some operator diagonal in momentum space, typically determined by the un- 
perturbed Hamiltonian, and Q is diagonal in coordinate space. The functional integral is 
taken with respect to some probability measure d/i(Q) and goes over the matrix elements 
of Q. [■ ]~y denotes the x, y-entry of the matrix [P + Q}" 1 . Our starting point is always a 
model in finite volume and with positive lattice spacing in which case the operator P + Q 
and the functional integral in (1.1) becomes huge- but finite-dimensional. In the end we 
take the infinite volume limit and, if wanted, the continuum limit. 

Our treatment is based on the following identity which is obtained by repeated ap- 
plication of the Feshbach formula (Lemma 3.2 below). It is proven in Theorem 3.3. Let 
B = (B kp ) ktPeM £ C NxN , M some index set, \M \ = N and let 



G(k):=[B-\ k (1.2) 



Then one has 



G(k) = (1.3) 

Bkk + ^( — l) r+1 E B kp2 Gk(p2)B P2P3 ■ ■ ■ B Pr _ lPr Gk P2 -p r _ 1 (pr)B Pr k 

r=2 p 2 ---pr£M\{k} 

where G kl ... kj (p) = [(B st ) s ^ M \ {kl ... k]} ]J is the p,p entry of the inverse of the matrix which 
is obtained from B by deleting the rows and columns labelled by ki, • • • , kj. In Section 2 
we apply this formula to a matrix of the form B = self adjoint + ie Id, which, for e ^ 0, 
has the property that all submatrices (B st ) Stt€ M\{k 1 --k j } are invertible. 

There is also a formula for the off-diagonal inverse matrix elements. It reads 

N 

[B- l ] kp = -G(k)B kp G k (p) + ^(-l)^ 1 E G(k)B kt3 G k (t 3 )B t3t4 • • • B trP G k ,,.. tr (p) 

r=3 t 3 ---t r eM\{k,p} 

(1.4) 

These formulae also hold in the case where the matrix B has a block structure B kp = 
(Bka,pr) where, say, a, r e {T, 1} are some spin variables. In that case the B kp are small 
matrices, the G kl ... kj (p) are matrices of the same size and the 1/- in (1.3) means inversion 
of matrices, see Theorem 3.3 below. 
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The paper is organized as follows. In the next section we demonstrate the method by 
applying it to the averaged Green function of the Anderson model. The Schwinger- Dyson 
equation for that model reads G^ 1 = Gq 1 + S(Go) where S(G ) is the sum of all two-legged 
one-particle irreducible diagrams. Application of (1.3) leads to an integral equation G~ L = 
Gq 1 + cr(G) where <r(G) is the sum of all two-legged graphs without two-legged subgraphs. 
The latter equation has two advantages. First, X is the sum of one-particle irreducible 
diagrams, but these diagrams may very well have two-legged subdiagrams and usually 
these are the diagrams which produce anomalously large contributions. And second, the 
propagator for a(G) is the interacting two point function G, which, for the Anderson 
model, is more regular than the free two point function G which is the propagator for 
the diagrams contributing to £(G ). More precisely, the series for <r(G) can be expected 
to be asymptotic, that is, its lowest order contributions are a good approximation if the 
coupling is small, but, usually, the series for T,(Gq) is not asymptotic. 

For the many-electron system and for the <p 4 model repeated application of (1.3,4) 
amounts to a resummation of two- and four-legged subgraphs. This is discussed in section 
4. In section 5 we discuss how our method is related to the integral equations which can 
be found in the literature. The proof of the inversion formula is given in section 3. 



2 Application to the Anderson Model 

Let coordinate space be a lattice of finite volume with periodic boundary conditions, lattice 
spacing \jM and volume [0,L] d : 

r = {x = JL(m, • • • , n d ) I o < m < ml - 1} = (±z) d /(LZ) d (2.1) 

Momentum space is given by 

M := T tt = {k = |(mv ,m d ) | < m, < ML - l} = (^Z) d /(27tMZ) d (2.2) 
We consider the averaged Green function of the Anderson model given by 

(G)(x,x) := J l-A-z + \V}^ xl dP(V) (2.3) 

where the random potential is Gaussian, 

dP(V)= ne-%. (2.4) 

Here z = E + is and A is the discrete Laplacian, 

d 

[-A - Z + W] x/X , = -M 2 ^ {$x',x+ei/M + bxi,x-ei/M ~ ^x',x) ~ Z 5 X)X > + \V X 8 X , X > (2.5) 



By taking the Fourier transform, one has 

(G)(x,x') = jgiid E e ik( - x '~ x \G)(k) (2.6) 



keM 



(G)(k) = j \a k 5 Kp + ^ d v k y l dP{v) (2.7) 

J RN <i L Viv j k k 



where N d = (ML) d = \T\ = \M\, dP(v) is given by (2.10) or (2.11) below, depending on 
whether N d is even or odd, and 



ak 



4M 2 ^sin 2 [4] -E-ie (2.8) 



The rigorous control of (G)(k) for small disorder A and energies inside the spectrum of 
the unperturbed Laplacian, E G [0,4M 2 ], in which case a k has a root if e — > 0, is still an 
open problem [AG,K,MPR,P,W]. It is expected that lim £ \ h m L->oo(G')(&;) = l/(a k — a k ) 
where Imo" = 0(A 2 ). 

The integration variables v q in (2.7) are given by the discrete Fourier transform of 
the V x . In particular, observe that, if F denotes the unitary matrix of discrete Fourier 
transform, the variables 

= (FV) 9 = VW E e ~ iq * V * = (f ) 1 TF E e ~"^ = (t) 1 ^ ( 2 - 9 ) 

would not have a limit if V x would be deterministic and cutoffs are removed, since the V q are 
the quantities which have a limit in that case. But since the V x are integration variables, 
we choose a unitary transform to keep the integration measure invariant. Observe also 
that v q is complex, v q = u q + iw q . Since V x is real, u_ q = u q and w_ q = —w q . In order to 
transform dP(V) to momentum space, we have to choose a set M + C M such that either 
q G A4 + or — q G M. + . If iV is odd, the only momentum with q = — q or w q = is q = 0. 
In that case <iP(l / ) becomes 



For even N we get 



dP(v) =e~ * n e -^ + ^ j =£2 (2.10) 
V2?r g e.M+ 71 



ffP(u) = c -^"n+<) rf "o^n n e -(ug+wg) j"gg!a (2.11) 

where go — is the unique nonzero momentum for which = 2irM(l, ■ ■ ■ , 1) — ^m. 
Now we apply the inversion formula (1.3) to the inverse matrix element in 



(£>(*■)=/ Ja^+^-J^PM (2.7) 
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We start with the 'two loop approximation', which we define by retaining only the r = 2 
term in the denominator of the right hand side of (1.3), 



-Dfcfc - 2^ pG M\{k} &kp<^k{P)&pk 



Thus, let 



-i 



G(k) := ya k 5 Kp + ^=v k _ p ^ ^ ^ = G(k; v, z) (2.13) 

In the infinite volume limit the spacing 2n / L of the momenta becomes arbitrary small. 
Hence, in computing an inverse matrix element, it should not matter whether a single 
column and row labelled by some momentum t is absent or not. In other words, in the 
infinite volume limit one should have 

G t (p) = G(p) for L -> oo (2.14) 

and similarly G tl ... tj (p) = G(p) as long as j is independent of the volume. We remark 
however that if the matrix has a block structure, say B = (B ka!PT ) with a, r e {f, 1} some 
spin variables, this structure has to be respected. That is, for a given momentum k all 
rows and columns labelled by k\, k[ have to be eliminated, since otherwise (2.14) may 
not be true. 

Thus the two loop approximation gives 

G(k) = r ^ (2.15) 

ak + 7n^ v ° ~ W £p#* Vk ~p G ^> v P~ k 

For large L, we can disregard the -j^ v o term. Introducing a k = cr k (v, z) according to 

G(k) =: — L_ , (2.16) 
Ofc — <J k 

we get 



p^k ^ ^ P ^ ^ 



and arrive at 



(G)(k) = j dP(v) (2.18) 

J n A 2 \^ W-p\ 

jV<2 

Now consider the infinite volume limit L — > oo or iV = ML — > oo. By the central 
limit theorem of probability (|f g | 2 — (|f g | 2 )) is, as a sum of independent random 

variables, normal distributed. Note that only a prefactor of is required for that 
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property. In particular, if F is some bounded function independent of N, sums which 
come with a prefactor of 1/N d like J2 q c q\ v g\ 2 can be substituted by their expectation 
value, 

Jim f F(^^c fc K| 2 )rfP(t;) = FQim ^^c fe (K| 2 )) (2.19) 
Therefore, in the two loop approximation, one obtains in the infinite volume limit 

(G)(k) = =: (2.20) 

„ * y (K- P i 2 ) ak ~ {(Tk) 

u k m 2^ _ A 2 v (K-t\ 2 ) 

where the quantity (cr^) satisfies the integral equation 

(\Vk-p\ 2 ) L^oo A 2 /" d d p (\v k - p \ 2 ) 



(o"fc) - J^d 



^ a p - (dp) Md J [Qt 2nM]d (2vr) d a p - (a p ) 



For a Gaussian distribution (|f g | 2 ) = 1 for all g such that (a k ) = (a) becomes independent 
of k. Thus we end up with 

^ W = A^^d -2{k\ T-i ■ 7T (2-22) 

where (a) is a solution of 

d d p 1 
'[0,**]' (2^F 4M2 Eti sm 2 [&] - * - (a) 

= X 2 f - (2 23) 

7 [0i2w]d (2vr) d 4M 2 Eti s in 2 [f ] - ^ - (<j) ' 

This equation is of course well known and one deduces from it that it generates a small 
imaginary part Imcx = 0(\ 2 ) if the energy E is within the spectrum of —A. 

We now add the higher loop terms (the terms for r > 2 in the denominator of (1.3)) 
to our discussion and give an interpretation in terms of Feynman diagrams. To make 
the volume factors more explicit, asume that the lattice spacing in coordinate space is 
1/M = 1 such that N — L. 

For the Anderson model, Feynman graphs may be obtained by brutally expanding 

/OO p 
[a k Sk,p + ^v k . p ] "J dP = J2 (C[VC] r ) kk dP (2.24) 

r=0 J 

oo „ 
= VW a k a P2 l -a Vr a k I V k-p2 V pa-P3 ' ' ' V Pr~k dP 



Figure 1: A string of particle lines with unpaired squiggles (dashed lines) 
For a given r, this may be represented as in figure 1 (c& := l/a^). 

The integral over the v gives a sum of (r — 1)!! terms where each term is a product of r/2 
Kroenecker-delta's, the terms for odd r vanish. If this is substituted in (2.24), the number 
of independent momentum sums is cut down to r/2 and each of the (r — 1)!! terms may 
be represented by a diagram 




Figure 2: Lowest order diagrams contributing to (2.24) for r = 2 and r = 4 

where, say, the value of the third diagram is given by Y\ n „ . For 

short: 

(G)(k) = sum of all two legged diagrams . (2.25) 

Since the value of a diagram depends on its subgraph structure, one distinguishes, in the 
easiest case, two types of diagrams. Diagrams with or without two-legged subdiagrams. 
Those diagrams with two-legged subgraphs usually produce anomalously large contribu- 
tions. They are devided further into the one-particle irreducible ones and the reducible 
ones. Thereby a diagram is called one-particle reducible if it can be made disconnected 
by cutting one solid or 'particle' line (no squiggle or dashed line), see also figure 3. 



(a) reducible (b) irrdeducible, no two-legged (c) irreducible, with two-legged 

subgraphs subgraphs 

Figure 3: Reducible and irreducible diagrams 

The reason for introducing reducible and irreducible diagrams is that the reducible ones 
can be easily resummed by writing down the Schwinger- Dyson equation which states that 
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if the self energy S fc is defined through 

<«>« - ^d^y < 2 - 26 > 

then Efc(Go) is the sum of all amputated (no l/o^'s at the ends) one particle irreducible 
diagrams. Here we wrote S fc (G ) to indicate that the factors ('propagators') assigned to 
the solid lines of the diagrams contributing to S fc are given by the free two point function 
Gq(p) = However, the diagrams contributing to X^Gq) still contain anomalously large 
contributions, namely irreducible graphs which contain two-legged subgraphs like diagram 
(c) in figure 3. 

In the following we show, using the inversion formula (1.3) including all higher loop 
terms, that all graphs with two-legged subgraphs can be eliminated or resummed by writing 
down the following integral equation for (G): 

(G)(k) = - -L_ (2.27) 

a k - Vk{{G)) 



where 



<Jk{{G)) is the sum of all amputated two-legged diagrams which do not contain 
any two-legged subdiagrams, but now with propagators (G)(k) = a J_ a ^ instead 
o 



Of Gn = f 



which may be formalized as in (2.35) below. The advantage of this formula is that the 
series for (((?)) can be expected to be asymptotic, that is, its lowest order contributions 
are a good approximation if the coupling is small, but, usually, the series for Efc(Go) is 
not asymptotic. Thus, in order to rigorously controll (G)(k), one has to define a suitable 
space of propagators, to estimate the sum of all amputated two-legged graphs without 
two-legged subgraphs on that space and then finally to show that the equation (2.27) has 
a solution on this space. We intend to address this problem in another paper. 

We now show (2.27) for the Anderson model. For fixed v one has 

G(k,v) = (2.28) 

a k - crk{v) 

where 

Vk(v) = J2(- 1 TY1 V!T G k(P2)--- G kp2---Pr-i(Pr)vk-p 2 ---v Pr -k (2.29) 

We cutoff the r-sum in (2.29) at some arbitrary but fixed order I < L d where £ is choosen 
to be independent of the volume. Furthermore we substitute Gk P2 ... Pj (p) by G{p). Thus 

(G)(k) = ( } ) (2.30) 

\a k - Lr=2 W)/ 
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where 

°k(v) = (-l) r jbr E G(p 2 )---G(p r )v k _ P2 ...v Pr _ k (2.31) 

Pi^Pj, Pi^k 

Consider first two strings s^ 1 , s r k 2 where 

s l( v ) = E c *P2-pt Ufe -P2 ' ' ' *V- fe ( 2 - 32 ) 

Pii'Pj, Pi^k 

and the c£ p2 ... Pr are some numbers. Then in the infinite volume limit 

<« 2 > = (4 2 ) (2-33) 

because all pairings which connect the two strings have an extra volume factor 1/L d . 
Namely, if the two strings are disconnected, there are {r\ +r 2 )/2 loops and a volume factor 

j (V H-T2 ) 

of 1/y L d giving (rx + r 2 )/2 Riemannian sums. If the two strings are connected, there 
are only (r x + r 2 — 2)/2 loops leaving an extra factor of 1/L d . By the same argument one 
has in the infinite volume limit 



(KT 1 ■ ■ ■ = (4T 1 ■ ■ ■ (sl m ) nm (2.34) 



which results in 



(G)(k) = (2.35) 

ak ~Yl VW Yl ( G )^) ' ' ' ( G )(Pr) (Vk-p 2 ■ --Vpr-k) 
r=2 P2---Pr 

Pi^Pj,Pi^k 

The condition p 2 , ■ ■ ■ , p r ^ k and Pi ^ pj means exactly that two-legged subgraphs are 
forbidden. Namely, for a two-legged subdiagram as in (c) in figure 3, the incomming and 
outgoing momenta p, p' (to which are assigned propagators (G)(p), (G)(p')) must be equal 
which is forbidden by the condidtion pi ^ pj in (2.35). 

However, we cannot take the limit £ — > oo in (2.35) since the series in the denominator 
of (2.35) is only an asymptotic one. To see this a bit more clearly suppose for the moment 
that there were no restrictions on the momentum sums. Then, if V = (^=gVk- p )k, P and 
(G) — ((G) (k) 8 k ,p)k,p, 



P2—Pr 



£ (G)(p 2 ) ■ ■ ■ (G)(p r ) (v k . P2 ■ ■ -v Pr - k ) = ((VmVY- 1 )^) (2.36) 



and for £ — > oo we would get 



{G)(k) ^-((^k) «*-((^v) ft ) (2 ' 37) 
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That is, the factorials produced by the number of diagrams in the denominator of (2.35) 
are basically the same as those in the expansion 

/ = £_^(2r + l)!! + R e+1 (\) (2.38) 

where the remainder satisfies the bound |i^ + i(A)| < t\ const^ \ e . 

We close this section with two further remarks. So far the computations were done 
in momentum space. One may wonder what one gets if the inversion formula (1.3) is 
applied to [—A + z + AV] -1 in coordinate space. Whereas a geometric series expansion of 
[— A + z + AV]^ 1 gives a representation in terms of the simple random walk, application 
of (1.3) results in a representation in terms of the self avoiding walk: 

[-A + * + AV%-i= V ^tV^r'^ 1 (2-39) 

■7 self avoiding 

where T is the lattice in coordinate space. Namely, if \x\ > 1, the inversion formula (1.4) 
for the off-diagonal elements gives 

L d 

[-A + \V\£ = ^T(-ir +1 J2 G(0)G (x 3 ) ■ ■■G 0xa ... Xr (x) (-A) 0X3 ■ • • (-A)^ 

r=3 x 3 ---x r er\{o,x} 

= det [(-A + W) y , y >er\{o}] det [(-A + \V)y, y 'er\{o,x 3 -x r ,x}] 

~ U * 2 =o„3,..™ +1 =* e r ^t [(-A + \V) y , y , er ] " det [(-A + \V) y , y , er \ {0 , x ,.. Xr} ] 

which coincides with (2.39). 

Finally we remark that, while the argument following (2.32) leads to a factorization 
property for on-diagonal elements in momentum space, {G{k) G{p)) = (G(k)) (G(p)), there 
is no such property for products of off-diagonal elements which appear in a quantity like 

A (<?) = h Yl \ \ a kSk, P + -^v k _ p ] \a k 5k, P + -^nv k ^ p ] ) (2.40) 

' \L V L J k,p L vi J fc— q t p— q 

k,p 

which is the Fourier transform of (| [— A + z + Al^]"* | 2 ). (Each off-diagonal inverse matrix 
element is proportional to 1/VI/, therefore the prefactor of 1/L d in (2.40) is correct.) 
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3 Proof of the Inversion Formula 



Lemma 3.1: Let B G C kxn , C G C nxk and let Id k denote the identity in C kxk . Then: 

(i) Id k — BC invertible Id n — CB invertible. 

(ii) If the left or the right hand side of (i) fullfilled, then C Id * r = Id l _ CB C. 



Proof: Let 

(' ^ 

B = j 

V - h 



C = ci • • • 4 



where the bj are n-component row vectors and the Cj are n-component column vectors. 

— » 

Let x G Kern(/<i — C5). Then x = CBx = Ylj^jCj if we define Aj := (bj,x). Let 
A = (Aj)i<j<fc- Then [(id - BC)X]i = \- J2j(hcj)Xj = {b u x) - Y.jih.c^Xj = since 

— * — * 

x = J2j thus ^ e Kern( Id — .BC). On the other hand, if some A G Kern(7d — BC), 
then x := ^jXj&j G Kern(/<i — C5) which proves (i). Part (ii) then follows from C = 
id n -CB ( Id n ~ CB)C = Idn \ CB C{Id k - BC) M 



The inversion formula (1.3,4) is obtained by iterative application of the next lemma, 
which states the Feshbach formula for finite dimensional matrices. For a more general 
version one may look in [BFS], Theorem 2.1. 



Lemma 3.2: Let h = (£ G C nxn where A G C kxk , D G c(n-k)x(n-k) are inverti bl e 

and B G C kx(n ~ k) , C G C (n - k)xk . Then 

h invertible A-BD~ X C invertible D-CA~ l B invertible (3.1) 

and if one of the conditions in (3.1) is fullfilled, one has h^ 1 = ^) where 

E = A-BD-^C ' H = D-CA- 1 B ' (^- 2 ) 

F = -EBD- 1 = -A~ X BH , G = -HCA^ 1 = -D~ 1 CE . (3.3) 



Proof: We have, using Lemma 3.1 in the second line, 

A-BD- l C inv. Id k - A^BD^C inv. 

Idn-k- D^CA^B inv. 
D-CA^B inv. 
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Furthermore, again by Lemma 3.1, 

D 1( ^ Id-A-^BD-^C = Id-D-WA-^B ® ^ = D-CA^B ^ = 

and 

A — 1 D 1 1 4-lp 1 JD jp T3 

71 D Id-D-^CA-^B ~ Id-A-^BD-^C 71 D ~ A-BD^C D ~ ^ D 

which proves the last equalities in (3.3), HCA' 1 = D~ X CE and EBD~ l = A^BH. Using 
these equations and the definition of E, E, G and H one computes 

(A B\ IE F\ _ IE F\ (A B\ _ (Id \ 
KC D) Kg h) ~ KG h) Kc d) ~ K Id) 

It remains to show that the invertibility of h implies the invertibility of A — BD~ 1 C. To 
this end let P = ( Id J,P= (° Id ) such that A — BD~ 1 C = PhP — PhP(PhPy 1 PhP . 
Then 

(A-BD- 1 C)PhT 1 P = PhPh' 1 P - PhP(PhP)- 1 PhPh' 1 P 

= Ph(l - P)/i _1 P - PhP(PhPy 1 Ph(l - P> _1 P 
= P - PhPh- 1 P + PhPh- 1 P = P 

and similarly Ph' 1 P(A - BD~ l C) = P which proves the invertibility of A - BD~ X C ■ 



Theorem 3.3: Let B G C nNxnN be given by B — (B kp ) kpe _\4, At some index set, 
| At | = N, and B kp = (B kapT )^ TeI G C nxn where I is another index set, |/| = n. Suppose 
that B and, for any M C At, the submatrix (B kp ) kjPe ^f is invertible. For k G At let 

G(fc):= (3.4) 

and, if A/ C At, k£N, 

G N {k) := [{(fl*).,*^}" 1 ]** G C" x " (3.5) 

Then one has 

(i) The on-diagonal block matrices of B~ l are given by 

G(k) = (3.6) 

B kk — y^(~l) r y^ B kp2 Gk{p2)B P2P3 ■ ■ ■ B Pr _ lPr Gkp 2 -p r _ 1 {pr)Bp r k 

r=2 p 2 ---PreM\{k} 

where 1/ • is inversion ofnxn matrices. 
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(ii) Let k,p E M, k ^ p. Then the off-diagonal block matrices of B 1 can be expressed 
in terms of the GV(s) and the B st , 



N 



[B-% p = -G(k)B kp G k (p) - ^T(-l) r G(k)B kt3 G k (t 3 )B t3U ■ ■ ■ B Up G kt3 ... tr (p) 



r=3 t 3 —t r eM\{k, P } 



(3.7) 



Proof: Let k be fixed and let p,p' E M. \ {k} below label columns and rows. By Lemma 
3.2 we have 



/ G(k) 
\ 

where 

G{k) = E 



\ 



( B kk — B kv — \ ( E — F — \ 

G H 



'kk t^kp 
B p ' k 



Bp/p 



\ 



J 



\ 



J 



Bkk~ ^2 Bk P [{{Bp'p)p', p eM\{k}} ^^Bp'k 

p,p'^k 

1 

Bkk — S ^B kp G k (p)Bp k — ^2 Bkp \_{{Bp<p)p',peM\{k}} 1 ] pp , Bp> k 

Py^k P.p'^k 

p¥=p' 



(3.8) 



and 



F kp =[B x ] = -G(k)^2B kt [{(B p/p ) p ' !peM \ {k} } : ] 



t P 



(3.9) 



t^k 



-G(k)B kp G k (p) - G(k) ^2 B kt [{(B p >p)p' tPe M\{k}} x ] 



Apply Lemma 3.2 now to the matrix {(-B p 'p)p>e.M\{fc}} 1 an d proceed by induction to 
obtain after £ steps 



G(k) 



B k k - y^(~l) r y^ B kp2 G k (p 2 )B P2P3 ■ ■ ■ Bp r _ 1 p r G kp2 ...p r _ 1 (p r )Bp r k - Re+i 

r=2 P 2 ---pr£M\{k} 
Pi^Pj 



(3.10) 



Fkp 



-G(k)B kp G k (p) - ]T(-l) r J2 G(k)B kt3 G k (t 3 )B t3t4 ■ ■ ■ B trP G kt3 ... tr (p) - R e+l 

(3.11) 



r=3 t 3 ---t r eM\{k, P } 
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where 



R e+1 - B kP2 G k{P2) • • -Bvt-xvi [{( B p'p)p',peM\{k P2 -pe}} 1 ww 5 K+i'! 



(3.12) 



R t+i - {~ l Y ^2 G ( k ) B kti ■ ■ ■G kta ... te _ 1 (t e )B tete+1 [{{B p ' p ) p ' ;PeM \{kt 3 -t e }} l ] tl+lP 

t 3 -t e+1 eM\{k :P } 

(3.13) 

Since Rn+i = Rn+i = the theorem follows ■ 



4 Application to the Many- Electron System and to 
the (/-Model 



4.1 The Many-Electron System 

We consider the many-electron system in the grand canonical ensemble in finite volume 
[0, L] d and at some small but positive temperature T — 1/(5 > with attractive delta- 
interaction given by the Hamiltonian 

H = H - \H int = ^(|^ - n)a£ a a ka - ^ a kT a ct-k| a q-pl a pt I 4 - 1 ) 

kcr kpq 

Our normalization conventions concerning the volume factors are such that the canonical 
anticommutation relations read {a ku) a^ T } = L d 5^-p5 u>T . The momentum sums range 
over some subset of say M = {k G | \e^\ < l}, e k = k 2 /2m — /i, and 

q G {k - p |k,p G M}. 

We are interested in the momentum distribution 

(a+ a a kCT ) = Tr[e^ H a+ a a^}/Tr e' pH (4.2) 

and in the expectation value of the energy 

(tf int > = £A(q) (4.3) 
q 
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where 

A(q) = Tr [e"^a+ aJ_ ki o q _ pi a pT ]/Tr e"^ (4.4) 

k,p 

By writing down the perturbation series for the partition function, rewriting it as a 
Grassmann integral 



Tre -/3(H -XH int ) 



(4.5) 



due = n^e-^^^^n^, 



performing a Hubbard-Stratonovich transformation (0 g = u q + ii> ? , d(p q d(p q := du q dv q ) 

e -E,«A _ I e iY. q {a q 4> q +b q ^ q ) e -Y. q \<t>q\ 2 j^^i± (4.6) 



with 



a„ = 



A? 



(4.7) 



and then integrating out the ■*/>,-?/> variables, one arrives at the following representation 



which is the starting point for our analysis (for more details, see [FKT] or [LI]): 

Jd(a^ a a ka ) = ± ^ (^kkoa^koa) 

feoef (2Z+i) 

where, abbreviating A; = (k, fc ), k = afc = iko — e^, g = A^, 



(4.8) 



dP(<j>) 



(4.9) 



ta,t<r 



and dP((f)) is the normalized measure 



dP(0) = 1 det 



— <P P -k 



k@k,p 



Furthermore 



where 



(4.10) 



(4.11) 



A(?) 



fc,p 



= (10, 



(4.12) 
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and the expectation in the last line is integration with respect to dP((j>). The expecta- 
tion on the ip variables {^ka^ka) = f i>koi>ka ^, P , q M%-u%- P ^ ^ c is Grassmann 
integration, but these representations are not used in the following. The matrix and the 
integral in (4.9) become finite dimensional if we choose some cutoff on the k variables 
which is removed in the end. The set M. for the spatial momenta is already finite since 
we have chosen a fixed UV-cuttoff |ek| = |k 2 /2m — fj,\ < 1 which will not be removed in 
the end since we are interested in the infrared properties at k 2 /2m = 

Our goal is to apply the inversion formula to the inverse matrix element in (4.9). Instead 
of writing the matrix in terms of four N x N blocks (a k 5 k ^ p ) k ^ p , {4> P -k)k,p-> {4>k-p)k,p and 
(a_fc5fc iP )fc iP where N is the number of the d+ 1-dimensional momenta k,p, we interchange 
rows and columns to rewrite it in terms of N blocks of size 2x2 (the matrix U in the 
next line interchanges the rows and columns): 

U 1 = B = (B kp ) k>p 
where the 2x2 blocks B kp are given by 

We want to compute the 2x2 matrix 

(G)(k) = J G{k)dP{4>) (4.14) 

where 

G(k) = [B-% k (4.15) 



U 



29 x 



12- 6 



We start again with the two loop approximation which retains only the r = 2 term in 
the denominator of (1.3). The result will be equation (4.20) below where the quantities 
(<7fc) and (|0o| 2 ) appearing in (4.20) have to satisfy the equations (4.21) and (4.24) which 
have to be solved in conjunction with (4.29). The solution to these equations is discussed 
below (4.30). 

We first derive (4.20). In the two loop approximation, 



G{k) 



Bkk — ^ B kp G k (p) B pk 

p+k 




l 9 sh 




4>p-k 



G k (p) 



4>k- P 



'} p -k 



(4.16) 
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where, substituting again G k {p) by G{p) in the infinite volume limit, 



4>p-K 



'k—p 



a 18- 

^|0o a, 



4>o 



+ S(p) 



'k—p 



'Jp-k 



(4.17) 



Anticipating the fact that the off-diagonal elements of (S) (k) will be zero (for 'zero external 
field'), we make the Ansatz 



E(fc) 



(4.18) 



and obtain 



-Y 

P¥=k 



(a>k + Vk)\<Pp~k\ 2 -^<pQ<pk- P <Pp-k\ , 



(a p +a p )(a p +* p )+±\M 2 \ <f) (f) k _ p <f) p _ k (d k + CT k )\(f) k _ p \ 



19) 



As for the Anderson model, we perform the functional integral by substituting the quan- 
tities \4> q \ 2 by their expectation values (l^l 2 )- Apparently this is less obvious in this case 
since dP((f>) is no longer Gaussian and the \(fi q \ 2 are no longer identically, independently 
distributed. We will comment on this after (4.37) below and at the end of the next sec- 
tion by reinterpreting this procedure as a resummation of diagrams. For now, we simply 
continue in this way. Then 



(G)(k) 



_i , a k + (a k ) (0 O > 

, Jfc )| 2+ A ( | 0o | 2> \ -i|(0 o ) a k + (a k ) 



where the quantity (ak) has to satisfy the equation 



p^k 



Since dP((f>) is not Gaussian, we do not know the expectations 
partial integration, we obtain 

p J 

Namely, 



(4.20) 
(4.21) 

W 2 ). However, by 
(4.22) 



g 2 > 



0, det [{B kp (<f>)} k!P \ e-^l^l 

% J ^ £ det 



i + 



ka,p'r' 



kcr,p"T" 
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Since 



we have 







Sk.-r 



det 



75 9Bka.pT TD 

J3ka,p'T' — ~ /■>/.■„ , 



>kcr,p"T n 



det [{B kp } Kp } = 



if r =| 

^[B-\ hp+ql ifr=t 



which results in (4.22). 

The inverse matrix element in (4.22) we compute again with (1.3,4) in the two loop 
approximation. Consider first the case q — 0. Then one gets 



1 + v*E / <t>oG{p) u dP{<t>) 

IE/ 



\a p +<r p \ 2 +±\<t> \ 2 



-^ Vo a p + a p 



dP{<f>) 



\a P +c P \ 2 +±\M 2 



dP{(j>) 



(4.23) 



Performing the functional integral by substitution of expectation values gives 



K+(^>l 2 +^(l4>ol 2 > 



or 



(l^ol 2 ) = 



i 

K./L^\a v +(a v )\ 2 - 



(4.24) 



K+K>l 2 +^(l^ol 2 > 



Before we discuss (4.24), we write down the equation for q ^ 0. In that case we use (1.4) 
to compute [-B -1 (0)] p -|-,p+gj. in the two loop approximation. We get 



[G(p)B PtP+q G(p + q)} 



11 



i i w_ x 

|a p +<T p | 2 +^|^o| 2 \a p+q +a p+q \ 2 + ±\<f> \ 2 V* 

+ ^P+qM-q + (a + a) p <f)o4> q } 

(a + a) p (a + a) p+q (f)_ q - ^4>l<P q 

ig_ {a+°')v(a+v) P +q4>q-^<i>o<t>-q 

(\a p +a p \ 2 +^\M 2 ) {K+o+T P+q \ 2 +^\M 2 ) 



(a + a) p (a + a) p+q (f> q - \^- q \ 
-^[(a + o-) P+g 0o0g + (a + a) p o 0_g] J 

(4-25) 



17 



which gives 



L {a+a)r{a+v)v+g<l>q--^4>o<l>-q 

9 (|ap+^P + ^|0ol 2 )(|a P+9 +^ +9 l 2 + ^l0ol 2 ) 



dP(</>) 



, A (ap + {av)){ap+q + {°P+q)){\<Pq\ 2 )-^{4>l<t>q4>-q) 

^ k (|a p + ( CTp >p + |(|^ p>)(|a p+9 +( ( T p+9 >P + A(|^|2}) 



(4.26) 



Although one may think that the expectation {<pQ<j) q <j)_ q ) vanishes for zero external field, 
this is not so. This can be seen again by partial integration: 

^>l<t> q <t>- q ) = iy^0^ det [{ Sfcp (0)} fcjp ] e-^\t>«\ 2 d<t> q d4> q 

= \ /^ 9 (flfcdet[{fl^)}*,p])e-S,l«.l a d^ 9 



D 9BkcT,pT p 

J-'k(T,p'T' d<f>- q - D ka,p"T" 



The above determinant is multiplied and devided by det [{-B fcp } fcp ] to give 



e-^^ 2 d(f) q d(f) q 



det 



-Dfccr,p'r' — oT^ £>ko,p"T" 



if r =t 



det [{-B fep } fei p] 



Computing the inverse matrix element again in the two loop approximation (4.25), we 
arrive at 

/ 72 1 I \ _ A V" 1 / ( a P+ a p){ap+q+ a P+q)4%<t>q4>-q-^<t>l<t>l<t>q<t>q 

\<P <P q <P- q ) K Z^\(|a p+CTp p + A| 0o | 2 )(| flp+9+CTp+9 |2 + A| (/)o | 2 ) 
P 

Abbreviating 

a P +(<r P ) f _ V^dM 2 ) 



9 P 



K+K>l 2 +^(l4>ol 2 >' 



fp 



K+(^>l 2 +^(l</>o| 2 > 



(4.27) 



this gives 



or 



A/ J2 a A \_ K Ylpfpfp+q K (l0o| 2 ) ^ ^ 



k 9p9p+q 

Substituting this in (4.26), we finally arrive at 

ni |2\ 1 ~~ K^p9p9p+q 



1 « XlpS'pS'p+gl ~ { K ^2 p fpfp+q) 



Jp o V ov+<l 
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(4.28) 



(4.29) 



where g p , f p are given by (4.27). Observe that, since dP(4>) is complex, also (\4> q \ 2 ) is in 
general complex. Only after summation over the q variables we obtain necessarily a real 
quantity which is given by (4.4,11). 

We now discuss the solutions to (4.24) and (4.29). We assume that the solution (a k ) 
of (4.21) is sufficiently small such that the BCS equation 

«E K+K>p+|Ap = 1 ( 4 - 30 ) 

p 

has a nonzero solution A ^ (in particular this excludes large corrections like (a p ) ~ , 
a < 1/2, which one may expect in the case of Luttinger liquid behaviour, for d = 1 one 
should make a seperate analysis), and make the Ansatz 

\(\<P \ 2 )=f3L d \A\ 2 + r ] (4.31) 

where 77 is independent of the volume. Then 

IV" I AT 1 

k |a p +(a p >|2+A(|0 o |2) k M^yFRAFTf 

V P 

« Z_> |a p +(a p )| 2 +|A|2 K (|a p + ( CTp >p + |A|2)2 ) 

= l-^ + 0((^) 2 ) (4.32) 

where we put ca = ^ X) p (| a +(0- )| 2 +|A| 2 ) 2 anc ^ use d the BCS equation (4.30) in the last line. 
Equation (4.24) becomes 

and has a solution 77 = A/(ca|A| 2 ). 

Now consider (\4> q \ 2 ) for small but nonzero q. In the limit q — > the denominator in 
(4.29) vanishes, or more precisely, is of order 0(1 /k) since 

1 _ £ _ £ 51 /p/p = |a p + (^)p 1 +A ( | 0o |2 ) = 0(l/«) 

P P P 

because of (4.32). If we assume the second derivatives of (a^ to be integrable (which should 
be the case for d — 3 and (\4> q \ 2 ) ~ 1/q 2 by virtue of (4.21)), then, since the denominator 
in (4.29) is an even function of q, the small q behaviour of (\4> q \ 2 ) is 1/q 2 . This agrees with 
the common expectations [FMRT,CFS,B]. Usually the behaviour of (\4> q \ 2 ) is infered from 
the second order Taylor expansion of the effective potential 



Kfr({0j) = ^|0 g | 2 -logdet 



r ig_ <P P -k 

ig <t>k- P c 



(4.33) 
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around its global minimum [L2] 



0f n = y/pP^Koe** (4.34) 
where the phase 9 of O is arbitrary. If one expands V e s up to second order in 



VA [ p 9 e^ forg^O 

one obtains [L2] 

v cS ({M) = Kni„ + 2A)(po- v^^) 2 + ^K + ^)p: 



(4.35) 



2 

9^0 



where for small q one has ct q ^ q ~ g 2 . Hence, if V e s is substituted by the right hand side 
of (4.36) one obtains (|0 g | 2 ) ~ 1/q 2 . 

For d — 3, this seems to be the right answer, but in lower dimensions one would expect 
an integrable singularity due to (4.21) and (4.3,4,11). In particular, we think it would be 
a very interesting problem to solve the integral equations (4.21,24,29) for d — 1 and to 
check the result for Luttinger liquid behaviour. A good warm up excercise would be to 
consider the + 1 dimensional problem, that is, we only have the ko,po, qo~ variables. In 
that case the 'bare BCS equation' 

I Yl peFW = 1 

poGf (2Z+1) 

still has a nonzero solution A for sufficiently small T = 1/(3 and the question would 
be whether the correction (a po ) is sufficiently big to destroy the gap. That is, does the 
'renormalized BCS equation' 



E 



\po + \ a r,,! -— A - 

P „e|(2Z+i) 



2\ 



(o~ P() ) being the solution to (4.21,24,29), still have a nonzero solution? We remark that, if 
the gap vanishes (for arbitrary dimension), then also the singularity of (|</> g | 2 ) disappears. 
Namely, if the gap equation has no solution, that is, if ^ J2 p ^ + ^ ^ 2 < oo, then (|0 O 
given by (4.24) is no longer macroscopic (for sufficiently small coupling) and ^(|0o| 2 
vanishes in the infinite volume limit. And the denominator in (4.29) becomes for q — > 

1 _ a i 

K / j \a p +(<T p )\ 2 
V 
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which would be nonzero (for sufficiently small coupling). 



Finally we argue why it is reasonable to substitute |0o| 2 by its expectation value while 
performing the functional integral. We may write the effective potential (4.33) as 



where 



Wo) 



Kff({0j)= Wo)+^2({0J) 



(4.37) 



= K 



\M 2 _ i 



g 



(4.38) 



and 



V 2 ({<M)=]T|^| 2 -logdet 



'k,p 



ig iox 



-1 



% al k ^,p ^,p 



$k,p 

ig 4>k- 



ig <Pp-k 



(4.39) 



If we ignore the o -dependence of V 2 , then the o -i n tegral 
/ F (i|0 o | 2 ) e-n^dfodfo JF( P 2 ) e-^osip) pd p 



J e- Vl( *Md(j) d(j) f e- KV vcs(p) p dp 

simply puts |0o| 2 at the global minimum of the (BCS) effective potential. 



(4.40) 



4.2 The (/-Model 

In this section we choose the y9 4 -model as a typical bosonic model to demonstrate our 
method. As in section 2, we start in finite volume [0, L] d on a lattice with lattice spacing 
1/M. The two point function is given by 



S(x,y) = {(p x (p v ) :-- 



f R Nd <p x <Py e 2 M 



9 l ^4 



tie E 1 ,„(-A+™ 2 ) I , s jJ I »Jj 



U x d(p x 



J RNd e~ E * ^ e -^E I , y (-A+m 2 ) I ,^^» Ux d(fa 



(4.41) 



where 



A + m 2 ) x ^ y = M d —M 2 ^2(S x ,y^ ei /M + 8x, y +ei/M ~ 26 x>y ) + m 2 S X:y 



1=1 



(4.42) 



21 



First we have to bring this into the form f[P + Q]~^ y dfi, P diagonal in momentum space, 
Q diagonal in coordinate space. This is done again by making a Hubbard Stratonovich 
transformation which in this case reads 



e 2 



IV a 2 



with 



x ^ 



(4.43) 



(4.44) 



The result is Gaussian in the y^-variables and the integral over these variables gives 

-l 

J x,y 



S(x,y) 



j^(-A + m 2 ) Xyy - -JjUu a .<S x J dP{u) 



(4.45) 



where 



dP(u) = ±det \j^(-A + m 2 ) x , y -^jL s u x 5 x , y \ \ ^*<\ldu 



(4.46) 



Since we have bosons, the determinant comes with a power of —1/2 which is the only 
difference compared to a fermionic system. In momentum space this reads (compare 
equations (2.7-11)) 

S{x-y) = i- d ^ k{x - y \G){k) 



where ( 7? = v q + iw q , 7 _ 9 = j q , dj q d% := dv q dw q ) 



(G)(k) 



a k $k, P ~ 7p7fc- P kk dP(l) 



and 



dP(-y) = ±det a k 5 kjP - -jfaik-p 'e ^dv IT e l7 ^d 7(? d 7(? 



-l7 ? l 2 



(4.47) 



(4.48) 



(4.49) 



and .M + again is a set such that either q E A4 + or— g G Furthermore 



a fc = 4M 2 £sin 2 [^] + 



m 



i=i 



(4.50) 



Equation (4.48) is our starting point. We apply (1.3) to the inverse matrix element in 
(4.48). In the two loop approximation one obtains (70 = vq e R) 



where 



kk 



a k -^ + h E P ^ G fc (p)| 7fc - P | 2 a* + 



9 2 V 



1k-p\ 



(4.51) 



(4.52) 
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which results in 



(G)(k) 



1 



a k + {(7k) 



(4.53) 



where {a k ) has to satisfy the equation 



(l7 fc - P | 2 ) 
°p + (^p) 



9 

2L d 



(17,-pI 2 ) g2 v-(l7 fc - P | 2 ) + ^ 



dp + (<Tp) 



9 



p^fc 



where the last line is due to 



(vo) 



i> det 



2 _1„2 



°p + (°p) 



a k h,p-^mlk- P 'e^dv Q IT e 179,2 d 7<? d 7(/ 



(4.54) 



/Id 



&det 



flfc5fc, P - yhik-p] 2 }e 5«8d Uo n e l79|2 rf7 g rf 7g 

= -5 £ (- A) / h 5 ^ - ^ 7fc - p ]p"p rfp(7) (4 - 55) 

p pp 

As for the Many-Electron system, we can derive an equation for (| 7(? | 2 ) by partial 
integration: 



I 2 
2 



(l7 9 | 2 ) = \ I 7<?7gdet a k 6 k ,p - -jhlk-p 'e " dv Ue W^cfry, 



= 1 + 



1 - 



Z j 7qg 9 lq 



det 



a k^k,p ^q^lk-p 



= 1 



2 v 7 ! 3 y 



TT-det 


dk^k,p — -^=slk-p 


det 


a k S k , P ~ -^ilk-p 





lq 



a k Sk, P - -jhlk-p dP(l) 

v u J p,p+q 



e "doolie 2 ' 79 ' 2 (i 7g (i 7g 
dP( 7 ) 

(4.56) 



Computing the inverse matrix element in (4.56) again in the two loop approximation, one 
arrives at 

(I7,| 2 ) = l-(|7 9 | 2 )^E 

(a p +(a p ))(a p+q + (cr p+q )) 

P 

or 

|2\ 



1 



1 _i_ si f dd p 

2 Jfo ^-Ti-n/fiti f9^7v2 



1 



[0,27rM] d (27r) d (a P + K»(a p+ , + ( CTj)+g )) 

which has to be solved in conjunction with 

( ai \- a 2f d d p (hk- P \ 2 ) + k 

\ u kl y J[ ,27rM] d (2^ a p +{a p ) 
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(4.57) 



(4.58) 



Introducing the rescaled quantities 

K) = M 2 S ^, (\ % \ 2 ) = \±, a k = M 2 e*, e k = £ sin 2 | + ^ (4.60) 



d 

,2 



= 1 



(4.57,58) read 



s k = M-V/ [o , 2 ^(0^f (4.61) 

\ = 5 ^ ( 4 - 62 ) 

1 ~r 2 J[0,27r]d (27r)d (e p +s p )(e p+9 +Sp +9 ) 

Unfortunately we cannot check this result with the rigorously proven triviality theorem 
since (a k ) and (|7 g | 2 ) only give information on the 2-point function S(x,y), (4.41), and 
on ^i^j^) 4 ) = J2qMl) where A(q) = (|7 g | 2 ) — 1. However, the triviality theo- 
rem [F,FFS] makes a statement on the connected 4-point function S4 >c (x 1 , x 2 , x 3 , #4) at 
noncoinciding arguments, namely that this function vanishes in the continuum limit in 
dimension d > 4. 

Before we include the higher loop terms of (1.3,4) and give an interpretation in terms 
of diagrams, we would like to comment shortly on a problem which was suggested to us 
by A. Sokal after a preprint of this paper was published on the web. It refers to the <j) 2 ip 2 - 
or 2 02" m odel. That is, we have two scalar bosonic fields on a lattice with unit lattice 
spacing with action 

2 

s(<Pu <t>2) = J2Yl <M^)(- A + ™ 2 )<j>i(x) + A 4>i{xf4>2{x) 2 

i=l x x 

The question is whether there is exponential decay (or a gap in momentum space) for the 
two point function G(x,y) = f <f)i( x )My) e ~ S{<l>) / I e ~ S{<t>) in the 

zero mass m — » limit. 
A computation with the above formalism in two loop approximation gives G(k) = 
where the gap a has to satisfy the equation a = A Jj ^ ^ d which gives 

0(A) ifd>3 
0(Alog[l/A]) if d= 2 
O(Al) ifd=l 

We now include the higher loop terms of (1.3,4) and give an interpretation in terms of 
diagrams. The exact equations for (G)(k) and (|7 g | 2 ) are 

(G)(k) = J [a k 5 k „-fa k ^Y*dPW = (^) (4.63) 

N d r 

° k = ~7h Vo + 'l2{^h) Yl G 'fc(P2)---G' fe p 2 ... Pr _ 1 (p r )7 fc _p 2 7p 2 „ P 3---7p r _fc 



r = 2 P2---Pr^k 
Pi^Pi 
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and 



(|7,P 



1 + 5& E / 7, a fc 5 fclP - ^7,- P + dP( 7 ) 

zvi^ ; Jp,p+q 



(4.64) 



P = 2 l+lEfel E (G , (P2)G P2 (p 3 )---G ?) ,.. Pr _ 1 (p r )G' p ,.. Pr (p 2 + g )x 

r=2 P2 , "Pi-^P2+9 



p 2 "-J>r^P2+9 

Pi^p-i 



7p2-P3 ' ' ' lPr-l-Prlpr-P2-ql'P2+q-P2 

For r > 2, we obtain terms (7^ • • -7^) whose connected contributions 



> 



>■ 



Figure 4: three and higher loop contributions 





are, in terms of the electron or </? 4 -lines, are at least six-legged. Since for the many- 
electron system and for the y9 4 -model (for d — 4) the relevant diagrams are two- and 
four-legged [FT,R], one may start with an approximation which ignores the connected 
r-loop contributions for r > 2. This is obtained by writing 



<7fci---7fe„> ~ <7fc 1 ---7fc„> 2 
where (the index '2' for 'retaining only two-loop contributions') 



(4.65) 



<7fei • • -7fc2„> 2 : = E (7fe CT i7fc CT2 ) • • • (7fc ff( 2„-i)7fc CT2 J = y 7fei • -'Ik** dP 2 {l) (4.66) 



pairings a 



if we define 



E 



\lq\ 



dP 2 ( 7 ) := e ^ q W> n 



t<|7,| 2 > 



(4.67) 



Substituting rfP by dP 2 in (4.63,64), we obtain a model which differs from the original 
model only by irrelevant contributions and for which we are able to write down a closed 
set of equations for the two-legged particle correlation function (G)(k) and the two-legged 
squiggle correlation function (|7 q | 2 ) by resumming all two-legged (particle and squiggle) 
subdiagrams. The exact equations of this model are 



(G)(k) 



a k S k , P - 7p7fc-p| kk dP 2 (n) 



\%\ 2 ) = 1 + 



:i\~l, d E/ 
V 



7 9 



a k$k, P - ^hlk-p 



p,p+q 



dP 2 {l) 



(4.68) 
(4.69) 



25 



and the resummation of the two-legged particle and squiggle subdiagrams is obtained 
by applying the inversion formula (1.3,4) to the inverse matrix elements in (4.68,69). A 
discussion similar to those of section 2 gives the following closed set of equations for the 
quantities (G)(k) and (I7J 2 ): 

<G>( ^^>- « 2 > = ttW (4 ' 70) 

where 

i 

= iE^)^ + E(7p) r E (G>(P2)---<G>(Pr)<7*- w 7 w - M "-7 Pr -*>2 

(4.71) 

= "iE(7fe) r E E (s q , P3+1 - P AG)(P2)---{G)(p r ){G)(p 2 + q )x 



p r=2 P2---Pr^k 

Pi^Pj 



r-1 



r=2 S=3 P2---Pr¥=P2+9 

Pi+Pi 



(7P2-P3 ' ' 'lPs-Ps+1 ' ' ' lpr-l-Prlpr-P2-q)2 ^ 

(4.72) 

In the last line we used that 7^ in (4.64) cannot contract to 7 P2 -p 3 or to 7 Pr _ P2 _ (? . If the 
expectations of the 7- fields on the right hand side of (4.71,72) are computed according 
to (4.66), one obtains the expansion into diagrams. The graphs contributing to (o^.) 
have exactly one string of particle lines, each line having (G) as propagator, and no 
particle loops (up to the tadpole diagram). Each squiggle corresponds to a factor (I7I 2 ). 
The diagrams contributing to (71) have exactly one particle loop, the propagators being 
again the interacting two point functions, (G) for the particle lines and (I7I 2 ) for the 
squiggles. In both cases there are no two-legged subdiagrams. However, although the 
equation (|7 ? | 2 ) = 1+ ^ ^ resums ladder or bubble diagrams (which is apparent from (4.57) 
or (4.26)) and more general four-legged particle subdiagrams if the terms for r > 4 in (4.72) 
are taken into account, the right hand side of (4.71,72) still contains diagrams with four- 
legged particle subdiagrams. Thus, the resummation of four-legged particle subdiagrams 
is only partially through the complete resummation of two-legged squiggle diagrams. Also 
observe that, in going from (4.68,69) to (4.70-72), we cut off the r-sum at some fixed order 
£ independent of the volume since we can only expect that the expansions are asymptotic 
ones, compare the discussion in section 2. 



26 



5 Concluding Remarks 



In the general case, without making the approximation (4.65), we expect the following 
picture for a generic quartic field theoretical model. Let G and Go be the interacting and 
free particle Greens function (one solid line goes in, one solid line goes out), and let D 
and .Do be the interacting and free interaction Greens function (one wavy line goes in, one 
wavy line goes out). Then we expect the following closed set of integral equations for G 
and D: . 

G= G 1 + a(G, D) ' D = D^ + niCD) ^ 
where a and n are the sum of all two legged diagrams without two legged (particle and 
wavy line) subdiagrams with propagators G and D (instead of G , D ). Thus (5.1) simply 
eliminates all two legged insertions by substituting them by the full propagators. For the 
Anderson model D = D = 1 and (5.1) reduces to (2.27,35). 

A variant of equations (5.1) has been derived on a more heuristic level in [CJT] and 
[LW]. Their integral equation (for example equation (40) of [LW]) reads 

G = — = l - (5.2) 

G^ + a{G,D Q ) 1 ; 

where a is the sum of all two legged diagrams without two legged particle insertions, 
with propagators G and D . Thus this equation does not resum two legged interaction 
subgraphs (one wavy line goes in, one wavy line goes out). However resummation of these 
diagrams corresponds to a partial resummation of four legged particle subgraphs (for 
example the second equation in (5.4) below resums bubble diagrams), and is necessary in 
order to get the right behaviour, in particular for the many-electron system. 

Another popular way of eliminating two legged subdiagrams (instead of using inte- 
gral equations) is the use of counterterms. The underlying combinatorial identity is the 
following one. Let 

Sty,$) = J dki, k G- \k)i> k + SUi>^) (5-3) 

be some action of a field theoretical model and let T(k) = T(G )(k) be the sum of all am- 
putated two legged particle diagrams without two legged particle subdiagrams, evaluated 
with the bare propagator Gq. Let 5S(xjj,xjj) = J dk^kT{k)ipk- Consider the model with 
action S — 8S. Then a p-point function of that model is given by the sum of all p-legged 
diagrams which do not contain any two legged particle subdiagrams, evaluated with the 
bare propagator Gq. In particular, by construction, the two point function of that model 
is exactly given by Go. Now, since the quadratic part of the model under consideration 
(given by the action S — 5S) should be given by the bare Greens function Gq 1 and the 
interacting Greens function is G, one is led to the equation G" 1 - T(G) = G 1 which 
coincides with (5.2). 
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Since the quantities a and ix in (5.1) are not explicitely given but merely are given by 
a sum of diagrams, we have to make an approximation in order to get a concrete set of 
integral equations which we can deal with. That is, we substitute a and tc by its lowest 
order contributions which leads to the system 

G{k) = G (k)-i + jdpD(p)G(k- P y D{q)= D (q)-i + jdpG(p)G(p + q) (5 ' 4) 

This corresponds to the use of (1.3,4) retaining only the r = 2 term. Thus we assume that 
the expansions for a and ir are asymptotic. A rigorous proof of that is of course a very 
difficult mathematical problem and this has not been adressed in this paper. Roughly 
one may expect this if each diagram contributing to a and n allows a const™ bound (no 
n\ and of course no divergent contributions). One may look in [FKLT] for an outline of 
proof for the many electron system with an anisotropic dispersion relation. In that case 
actually one obtains a series with a small positive radius of convergence instead of only an 
asymptotic one (because the model is fermionic), which simplifies the proof considerably. 

Finally we remark that the equations (5.4) can be found in the literature. Usually they 
are derived from the Schwinger-Dyson equations which is the following non closed set of 
two equations for the three unknown functions G, D and T, T being the vertex function 
(see, for example, [AGD]): 

G(k) = G (k) + G (k) J dpG(p)D(k-p)T(p,k-p)G(k) 

D(q) = D (q) + D (q) J dpG(p)G(p + q)T(p + q,-q) D(q) (5.5) 

The function T(p, q) corresponds to an off-diagonal inverse matrix element as it shows 
up for example in (4.22). Then application of (1.4) transforms (5.5) into (5.1). One 
may say that although the equations (5.4) are known, usually they are not really taken 
seriously. For our opinion this is due to two reasons. First of all these equations, being 
highly nonlinear, are not easy to solve. In particular, for models involving condensation 
phenomena like superconductivity or Bose-Einstein condensation, it seems to be apropriate 
to write them down in finite volume since some quantities may become macroscopic. And 
second, since they are usually derived from (5.5) by putting V equal to 1 (or actually 
-1, by the choice of signs in (5.5)), one may feel pretty suspicious about the validity of 
that approximation. The equations (5.1) tell us that this is a good approximation if the 
expansions for a and re are asymptotic. 

The applications of the method shown in this paper basically confirmed the common 
expectations for the particular models, thus one may say there are no really new results. 
However, we think it is fair to say that the computation of field theoretical correlation 
functions is an extremely difficult mathematical problem and therefore one should have 
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welcome everything which sheds some new light on these problems. We hope that we 
could convince the reader that the method presented in this paper definitely does this. 
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